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1 Introduction 


This is the semi-annual status report for NASA grant NAG2-513, “The Sym- 
bolic Computation and Automated Analysis of Trajectories.” This report 
covers the period July 1, 1990 through December 31, 1990. The research 
supported by this grant is broadly concerned with the symbolic computa- 
tion and mixed numeric-symbolic computation of trajectories of dynamical 
systems, especially control systems. The NASA Technical Officer for this 
grant is Dr. George Meyer, NASA Ames Research Center, Mail Stop 210-3, 
Moffet Field, California, 94035. In Section 2, we review the progress dur- 
ing the past six months. Section 3 contains bibliographic references for the 
articles related to this grant that were completed during this period. 

2 Review of Work During Report Period 

2.1 Symbolic computation of vector fields and analytic algo- 
rithms 

This section is the introduction to [1], with a few minor modifications. 

During the report period, we further developed the algorithms in [7] and 
[10] for rewriting expressions involving differential operators. The differen- 
tial operators that we have in mind arise in the local analysis of nonlinear 
control systems. During this period, we extended these algorithms in three 
different directions: 

1. We generalize the algorithms so that they apply to differential opera- 
tors on groups. This generalization is important for applications. For 
example, the nonlinear system describing a robotic joint or a satellite 
evolves on the group G — 50(3) of spatial rotations. The local study 
of such systems requires the computation of expressions consisting of 
differential operators on G . 

2. We develop the data structures and algorithms to compute symboli- 
cally the action of differential operators on functions. Again, this is 
crucial for applications. For example, it is crucial to our approach 
for deriving the conditions for a numerical algorithm to remain con- 
strained to a group. In other words, if z n +\ = T(x n ) is the update 
rule for a numerical algorithm evolving on a group G, we would like 
to choose T so that x n G G implies x n +i G G. 
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3 . We showed that the algebraic formalism uderlying computations of 
this type can be given the structure of a quantum group. 

For a further discussion of applications of these algorithms, see [ 4 ] and [6] 
and the references given there. 

Here is the 

Setup. 

1 . Let k denote either the real or complex numbers. 

2. Let G denote a finite dimensional Lie group over k y g denote its Lie 
algebra, and Yi, . . Ys a basis for g of left-invariant vector fields. 

3 . Let R = C°°(G) denote the algebra of smooth functions on G taking 
values in k. 

4 . Fix M derivations of R of the form 

jV 

F ; = E a i y *. (i) 

and let A denote the free associative algebra k<F\ f . . F\f> of dif- 
ferential operators generated by F\, . . ., Fa/, with coefficients from 
Jb. 

We are concerned with the following 

Problem. Given a differential operator p £ A and a function f £ R } 
substitute the Equations ( 1 ) and compute p * / using as few operations as 
possible. This problem is interesting since in many cases cancellations take 
place. 

Example 1. Let G = denote the abelian group, 



the (left invariant) coordinate vector fields, and Fi, F2, F3 three fixed vec- 
tor fields defined in terms of the via Equations ( 1 ). Then the naive 
substitution of (1) and simplification ofp* /, where 

p = F3F2F1 — F3F1F2 — F2F1F3 + F1F2F3 £ v 4 , / G R, 

yields 24 N 3 terms, while more specialized algorithms need only compute the 
6 jV 3 terms which don’t cancel. These types of examples are considered in 
[ 7 ] and [ 10 ]. 
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Example 2. Consider the local analysis of a nonlinear system of the form 


x(t) = F(z(t)), r(0) = *°eG, (2) 


where 


M 

F = Y.», F : 

i=i 


In practice, the uj are constants, functions of time, or perturbation parame- 
ters. The study of this system typically involves the computations of various 
series in the algebra A of differential operators. For example, the local flow 
of the system is determined by the Taylor series 


exp hF = 1 + hF + ^F 2 + ^F 3 + ...€ A[{h]]. 

An alternative to computing higher derivatives F* is to choose constants 
c t , c,j, i = 1, . . ., fc, j < i, so that the expression 

exp hc^Fk * ' *exp hc\Fi } 

where 

N 

Fi = 53 aM ( x °) y M € 0 

M=1 

N 

a M (exp(/ic 2 i A ) * * 0 )F m G ^ 

N 

F 3 = ^(exp(ftc3 2 F 2 ) ■ exp(/ic 3 iFi) * x°)y /i € g , 


is equal to exp hF to order fc. Notice that the left invariant vector fields Fj 
arise by "freezing the coefficients” of F at various points along its flow. 

Expanding these expressions around the common base point x° E G 
yields many terms, which must cancel in the end if the algorithm is going 
to approximate the flow of the underlying nonlinear system. The action 
of the differential operators Fj on the coefficient functions a£ must also be 
computed. Notice, that unlike Example 1, the u here do not commute. 

The computations in both examples are easily kept track of by using 
finite rooted trees, labeled with the symbols Fi, . . ., F\f . It turns out the 
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the vector space, with basis the set of such trees, has an algebraic structure 
B which is crucial to efficiently organizing the computation. The advantage 
of working with the trees B is that many terms which cancel in the end need 
not be computed. See [6] for an expository treatment of this idea. The key 
observation required for this work is that it is possible to define an action of 
the algebra B of finite rooted trees, labeled with Fi, . . Fm , on the ring of 
functions R which enjoys essentially all the properties of the familiar action 
of the algebra A of differential operators generated by Fi, . Fm on R. It 
turns out that B is a Hopf algebra, just as A is, and that both actions give 
R the structure of what is called an H - module algebra. It also turns out 
that it is possible to define an alternative structure on both A and B which 
gives these algebras the structure of a quantum group. 

2.2 Very Large Databases of Control Trajectories 

During the reporting period, we completed an initial design of the system 
architecture for software to analyze nonlinear control systems using database 
computing, as described in [3] and [5]. The system will support a variety of 
objects suitable for analyzing nonlinear systems, including 

1. trajectory segment 

2. state 

3. and control system. 

The software will consist of the following components: 

L a storage manager responsible for managing the storage, retrieval and 
update of the objects supported by the system; 

2. an object manager responsible for keeping track of the various objects 
' recognized by the system and their relations; 

3. a query manager responsible for parsing queries and generating the 
necessary code for the object manager so that the appropriate objects 
can be retrieved. 

We expect software providing a “proof of concept” for this approach to be 
completed sometime during the coming year. 
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3 Articles Written During the Reporting Period 

In this section we include the title pages of the articles resulting from the 
research supported by this grant during the reporting period. 


7 



Using Trees To Compute Approximate Solutions to 
Ordinary Differential Equations Exactly 

October, 1990 


Abstract 

In this paper, we review some recent work relating families of trees 
to symbolic algorithms for the exact computation of series which ap- 
proximate solutions of ordinary differential equations. It turns out that 
the vector space whose basis is the set of finite, rooted trees carries a 
natural multiplication related to the composition of differential opera- 
tors, making the space of trees an algebra. This algebraic structure can 
be exploited to yield a variety of algorithms for manipulating vector 
fields and the series and algebras they generate. 


Acknowledgments 

Robert Grossman was supported in part by NASA Grant NAG 2-513 and 
by NSF Grant DMS-8904740. 

Status 

Computer Algebra and Differential Equations , M. F. Singer, editor, Aca- 
demic Press, New York, 1991, in press. 
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Bialgebra deformations of certain universal 
enveloping algebras 

Robert Grossman, University of Illinois at Chicago 
Dave Radford, University of Illinois at Chicago 

October, 1990 


Abstract 

Let A denote a bialgebra over a field k and let A t = >l{[f]] denote 
the ring of formal power series with coefficients in A. Assume that A 
is a free algebra over k with a basis of primitives. We give a simple 
construction which makes A, a bialgebra deformation of A. Usually A t 
is neither commutative nor cocommutative. This construction yields 
deformations of bialgebras associated with families of trees. 


Acknowledgments 

Robert Grossman was supported in part by NASA Grant NAG 2-513 and 
by NSF Grant DMS-8904740. 

Status 

Submitted for publication. 
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Computations involving differential operators and 
their actions on functions 

Peter Crouch, Arizona State University 
Robert Grossman, University of Illinois at Chicago 
Richard G. Larson, University of Illinois at Chicago 

December, 1990 


Abstract 

In this paper, we further develop the algorithms described in the pa- 
per “The symbolic computation of derivations using labeled trees” R. 
Grossman and R. Larson, J. of Symbolic Computation , to appear, for 
rewriting expressions involving differential operators. The differential 
operators that we have in mind arise in the local analysis of nonlinear 
dynamical systems. In this work, we extend these algorithms in two 
different directions: We generalize the algorithms so that they apply to 
differential operators on groups and we develop the data structures and 
algorithms to compute symbolically the action of differential operators 
on functions. Both of these generalizations are needed for applications. 
This paper is preliminary: a final paper containing proofs and a further 
analysis of the algorithm will appear elsewhere. 


Acknowledgments 

Peter Crouch was supported in part by the NSF grant INT-8914643. Robert 
Grossman was supported in part by NASA Grant NAG 2-513 and by NSF 
Grant DMS-8904740. Richard Larson was supported by NSF Grant DMS- 
8904740. 

Status 

Proceedings of 1991 International Symposium on Symbolic and Algebraic 
Computaiion } ACM, 1991, in press. 
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Using Trees To Compute Approximate Solutions to 
Ordinary Differential Equations Exactly 


Robert Grossman" 
University of Illinois at Chicago 

October, 1990 


Abstract 

In this paper, we review some recent work relating families of trees 
to symbolic algorithms for the exact computation of series which ap- 
proximate solutions of ordinary differential equations. It turns out that 
the vector space whose basis is the set of finite, rooted trees carries a 
natural multiplication related to the composition of differential opera- 
tors, making the space of trees an algebra. This algebraic structure can 
be exploited to yield a variety of algorithms for manipulating vector 
fields and the series and algebras they generate. 

1 Introduction 

In this paper, we review some recent work relating families of trees to sym- 
bolic algorithms for the exact computation of series which approximate so- 
lutions of ordinary differential equations. It turns out that the vector space 
whose basis is the set of finite, rooted trees carries a natural multiplication 
related to the composition of differential operators, making the space of trees 
an algebra. This algebraic structure can be exploited to yield a variety of 
algorithms for manipulating vector fields and the series and algebras they 
generate. , 

In Section 3, we introduce and explore the algebraic structure of trees. 
Section 4 describes a simplification algorithm for the rewriting of symbolic 

‘Laboratory for Advanced Computing, Department of Mathematics, Statistics, and 
Computer Science, Mail Code 249, University of Illinois, Chicago IL 60680, gross- 
man &uicbert. eecs.uic.edu. This research was supported in part by the grants NASA 
N AG2-513 and NSF DMS-8904740. 
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expressions involving vector fields. Section 5 describes an algorithm for 
generating explicitly integrable flows associated with nilpotent Lie algebras. 
Section 6 exploits the relation between Taylor series and trees to study a 
class of intrinsic numerical integrators. We begin in Section 2 with some 
background. 

The results surveyed here are the work of a variety of mathematicians: 

1 especially want to mention the contributions of my collaborators Peter 
Crouch, Matthew Grayson and Richard Larson. The work on algebras of 
trees and its applications to symbolic computation is joint work with Richard 
Larson. All of the algorithms described here rest upon this foundation. The 
work on explicitly integrable flows and nilpotent Lie algebras is joint work 
with Matthew Grayson. The work on numerical algorithms evolving on 
groups is joint work with Peter Crouch. 

2 Background 

Consider a differential equation 

x(0=Ei(*(0) + “£i(*(O). *(0) = x°gr n , U) 

where E \ and Ei are vector fields and u is a parameter. In applications, 
u will b e either a small perturbation u = c, a control t — *> u(t), or simply 
the constant u = 1. Unless the vector fields E\ and E 2 are very special, no 
algorithm is known which will return the general solution to the system in 
closed form. Our objective is to find efficient algorithms to compute various 
approximate solutions of the differential equation exactly using symbolic 
computation. 

Although the impact of symbolic computation in this area is recent, the 
connection between the existence of closed form solutions and the approxi- 
mation of general solutions is a traditional theme, dating back to at least the 
the nineteenth century. One can distinguish two approaches. One, cham- 
pioned by Lie, is based upon algebra and geometry and concerns us here; 
the other, championed by Weirstrass and Poincare, is based upon complex 
function theory. * 

Consider a group of transformations acting on RA of the form 

^ j • Z n — / (xj , * . . , XjV , S j , . . . , S r ), C • ■ 1 

with the property that the group permutes the solutions of the nonlinear 
system (1). Lie asked the question [51] and [52], How can information about 
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the Iranformatton group be used to help integrate the differential equation? 
To answer this question. Lie introduced the tnfinitesimal generators of the 
group 

1 < £ < r 

k[ } e, ds k dz - - 

and showed that the A t satisfy 


r 

[A,, A,} = Y 

k—\ 

where [-, ] is the commutator, or Lie bracket , 

[Ai,A,] = AiAj -AjAt, 

and the c* are constants. For example, Lie showed that if there is a one 
parameter group of transformations permuting the solutions of a nonlinear 
system in the plane (x,t), then the integrating factor for the equation may 
be read off from the infinitesimal generator. 

Since Lie s time, this basic question has contributed to the development 
of a number of different fields: 

• The vector fields Aj generate a filtered Lie algebra, which is usually 
infinite dimensional, and is the infinitesimal version of the continu- 
ous pseudogroup of transformations generated by the <J>j. Prior work 
has focused on the geometry and structure theory of these algebras; 
important contributions have been made by Guillemin and Sternberg 
[37] and [36], and Hermann [52], [53], building upon the earlier work 
of Cartan, Ehresmann and Spencer. 

• Formal sums of iterated powers of vector fields, or Lie series, have 
been developed by Grobner [25], [26] and Knapp and Wanner [46], 
[47] into an operational calculus and used to approximate the solu- 
tions of differential equations. Lie transform methods have also been 
used in perturbation theory by Rand [59] and Meyer [54], in celestial 
mechanics by Deprit [18); and in particle physics by Dragt [20]. 

• Explicit series computations of solutions of differential equations have 
a number of interesting connections with combinatorics. Chen [13] 
makes uses of the shuffle product, Joyal [45], Labelle [49], Leroux [50], 
and Viennot [71] employ trees and species, while Rota, Krahaner and 
Odlyzko [62] exploit the umbral calculus. 
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• Ritt and Kolchin, building upon earlier work of Picard and Vessiot, 
developed the field of differential Galois theory. The goal is to obtain a 
theory describing the solvability of differential equations analogous to 
Galois' theory describing the solvability of algebraic equations. The 
survey by Singer [67] provides a good description of this field from 
the viewpoint of symbolic computation. Other relevant contributions 
include the beginnings of a differerential Grobner theory [3], [38], anal- 
ogous to the Grobner theory in commutative algebra; and a Hopf al- 
gebraic interpretation of Picard- Vessiot theory by Takeuchi [70]. 

• There is now a resurgence of interest in using symmetry groups to 
help integrate differential equations. This direction of research has 
been active in the Soviet Union for some time, but during the past 
decade there has been increased interest in the United States. Impor- 
tant contributions have been made by Olver [57], Schwarz [64], and 
Dluman and Cole [5]. Closely related is the study by Caviness [66] of 
conservation laws for differential equations. 

During the past several years, there has been increasing interest in sym- 
bolic computation and differential equations. Work has proceeded in a num- 
ber of directions, and is based upon both the Lie and the Weifstrass and 
Poincare traditions: 

• Zippel [73] is writing a modular symbolic computation system which 
supports the ability to call high quality numerical routines. Using such 
a system, he has shown how symbolic algorithms can be used to select 
appropriate numerical algorithms. 

• Guckenhemier [35] has produced the program kaos which allows the 
user to access a variety of algorithms to investigate a differential equa- 
tion from the point of view of modern dynamical systems. 

• There are a number of programs to compute symmetry groups of differ- 
ential equations, including one written by Char [12], and the programs 
SODE and SPDE written by Schwarz [64]. 

• Abelson and Sussman ar/d their group at MIT [1], [2] have combined 
techniques in artificial intelligence to produce software which automat- 
ically analyzes the qualitative features of a differential equation. 

• Wang [72] and Steinberg [68] have developed systems which use sym- 
bolic computation to produce high quality, optimized Fortran code to 
solve differential equations. 
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El t— > 


E2 I — > 


E3 I — > 


-'El 


'E2 


'E3 


Figure 1: The trees associated with the vector fields E } . 


• Della-Dora and Tournier [17] have used the fundamental ideas of Ramis 
[58] to produce a system to analyze linear ordinary differential equa- 
tions. They are now turning their attention to nonlinear systems. 

Point of view. The point of view taken here is to focus on the algorithmic 
aspects of the computation of the vector fields A j and their brackets and 
to use this information to develop appropriate algorithms which use exact 
symbolic techniques to approximately integrate the trajectories of the dif- 
ferential equation. As will become clear, there are a number of interesting 
points of contact between this approach and the approaches just described. 

In the following sections, we review data structures and algorithms for 
the symbolic computation of the flows of vector fields, and for the symbolic 
approximation of general flows by flows which can be studied symbolically. 

3 Vector fields and the algebra of Cayley trees 

In this section, we describe a data structure which is central to the algo- 
rithms we give for the symbolic computation of series which approximate 
the solutions of differential equations. The basic idea is to assign trees to 
vector fields as illustrated in Figure 1, and then to impose a multiplication 
on trees which is compatible with the composition of vector fields. 

Consider three vector fields 

E\ — a\D\ + • • • a^D^, E? = b\D\ + ■ ■ ■ bpj D\ , 


£3 = cjDi + ...csDn 

where D, = d/dn, and a,, b, and c, are smooth functions on R N Now 

E^ • E\ — 'y b-, ( Dj a, ) Dj + 'y bjd t Dj Dj 


and £3 ■ £2 • £1 * s equal to 

Ck{Dkbj){Dja,)Di + 'Y^Ckb J {DkD ] ai)D, + ^2ckb } (Dja,)DkDi 
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Figure 2: The trees associated with Equation 2. 


+ c^bjd, Dk Dj D, + ckb](Dkai) Dj D x + cj \c(Dkb } )a, D } D,. (2) 

Here the sum is for i,j,k = and hence involves 0(N 3 ) differentia- 

tions. It is convenient to keep track of the terms that arise in this way using 
labeled trees: we indicate in Figure 2 the trees that are associated with the 
six sums in this expression. 

An expression such as 

(£ 3 , (£ 2l £,]] = £ 3 £ 2 £i - £ 3 £i £2 - E 2 E l E 3 + E X E 2 E 3 (3) 

gives rise in this fashion to 24 trees corresponding to the 24 N 3 differenti- 
ations that a naive computation of this expression requires. On the other 
hand, 18 of the trees cancel, saving us from computing ISN 3 terms. We 
are left with 6 N 3 terms of the form (junir)D Ml . A careful examination of 
this correspondence between labeled trees and expressions involving the E } s 
shows that the composition of the vector fields Ej’s, viewed as first order 
differential operators, corresponds to a multiplication on trees. This multi- 
plication is illustrated in Figure 3. It turns out that this construction yields 
an algebra, which we call the ajgebra of Cayley trees. 

Here is a more precise description for the specialist, following [27] and 
[30]. Let k denote a field of characteristic 0. We say that a finite rooted 
tree is labeled with {Ei, . . . ,E\f) in case each node, except for the root, is 
assigned a element from this set. Let CT{E \ , ..., £w) denote the set of 
labeled trees and let fc{£T(£i, . Em)} denote the vector space whose 
basis consists of labeled trees in £T(£i, ..., Em) ■ Suppose that ti, 
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Figure 3: An example of multiplying two trees. 

*2 e CT( E \ , . . . , Em) are trees. Let S\ t . . . , s r be the children of the root 
of it. If has n + 1 nodes (counting the root), there are (n + l) r ways to 
attach the r subtrees of t\ which have s\ , . . . , s r as roots to the tree <2 by 
making each s t the child of some node of < 2 - The product * 1*2 is defined 
to be the sum of these (n + l) r trees. It can be shown that this product is 
associative, and that the trivial tree consisting only of the root is a right and 
left unit for this product. In [27], we define a comultiplication on k{CT(E u 
. , . , E \{ ) } and show that 

Theorem 3.1 The vector space k{CT(E\ t . . . , Em)} basts all equiva- 
lence classes of finite rooted trees is a cocommutative graded connected Hopf 
algebra. 

We call this algebra the algebra of Cayley trees generated by the set of 
labeled trees CT(E\ y . . . , Em) The relation between trees and differential 
operators goes back to Cayley [8], [9]. The coalgebra structure on this space 
is very similar to the coalgebra structure defined by Joni and Flota [44]. 
However, the coalgebra structure defined there was defined for individual 
combinatorial objects, rather than for a class of objects such as the family 
of rooted trees. Butcher [6] and [7] has also defined a multiplication on the 
vector space which is dual to the space of trees. This multiplication is closely 
related to the one just defined. 

4 Symbolic evaluation of vector field expressions 

When expressions involving vector fields, such as Lie brackets, are written 
out in coordinates, there is typically a lot of cancellation. Similar cancel- 
lation occurs in expressions involving Poisson brackets and when flows are 
concatenated, as in Campbell-Baker-Hausdorf expansions. In this section, 
we use the algebra of Cayley trees to exploit this cancellation in order to 
compute more efficiently formal expressions involving vector fields. 

This is the set up. Let R denote a ring of functions. In applications, R 
is usually either the ring of polynomial functions, rational functions, or C°° 
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functions. Fix several first order differentia) operators with coefficients from 

R 

N 

Ej = E a > D M. < 6 R, j = 1 A# (4) 

that are defined in terms of a basis of first order differential operators 

= - — , /i = 1, . . . , AT. 

Simplifying any expression in the E^’s using Equation (4) yields a differential 
operator, which we can view as an element of End R. 

We now define a homomorphism from the algebra of expressions in the 
Ej s to the algebra of trees. Indeed, the assignment in Figure 1 extends to 
an algebra homomorphism from the free associate algebra in the symbols E : 
to the algebra of Cayley trees k{CT(E\ , . . . , Em)} . It is straightforward to 
define a downward-pointing arrow so that the following diagram commutes: 

k<E Xt ...,E M > - k{CT{E lt ...,E M )} 

\ l (5) 

End R 


Algorithm 4.1 To rewrite expressions in the first order differential opera- 
tors Ej in terms of the basis 


d d 2 

dx^ 0 x h j 0 x 


P 1 1 M2 , - - = 1 , . > AT, 


compute the composition of the rightward and downward pointing arrows in 
the diagram above. 


In [28], [29] and [33], we show that the algorithm is much more efficient than 
naive substitution, which corresponds to computing the diagonal arrow di- 
rectly. In some common cases, the improvement in efficiency is exponential. 
We have implemented this and related algorithms in Maple, Mathematica 
and Snobol. 

We illustrate this algorithm by working the example of the last section 
following [31]: consider a higher order derivation of the form 


p — E 3 E 2 Ei — E 3 E\E 2 — E 2 E\E 3 + E\E 2 E 3 . 

Naive simplification requires computing 24 jV 3 terms of the form described 
in Table 1. The image of p in the algebra of Cayley trees contains 24 
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No. of terms 

Form of terms 

8A^ ' 

coeff. D M1 

12 N 3 

coeff. 

AN 3 

coeff. 


Table 1: Naive computation of the differential operator corresponding to p. 



Figure 4: The surviving labeled trees. 


trees, six for each of the four terms of p. For example, the six labeled trees 
corresponding to the first term are given in Figure 2. Eighteen of these 
trees cancel, leaving the six trees in Figure 4. The corresponding differential 
operator is equal to 

+ Y2 a 3 3 a-2 3 ( a^ 1 ) D U1 - ^2 a^ 3 a« ( D D n a? )D n , 


and contains 18 jV 3 fewer terms/ of the form indicated in Table 2 than does 
the naive computation of p. An example of the cancellation of labeled trees 
is given in Figure 5. 

To summarize: we have defined an algebraic structure on families of trees 
which mirrors the algebraic structure of formal expressions in the variables 
Ej t but which alleviates the need for computing intermediate expressions 
which cancel when the noncommuting Ej’s are expressed in terms of the 
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0 



Figure 5: The term E$E\E 2 contributes the first labeled tree and the term 
El E 2 £3 contributes the second, which cancel. 


No. of terms 

Form of terms 

2jV 2 3 

coeff. D mj 

\2N 3 

coeff. 

4N 3 

coeff. D ti ,D^D llt 


Table 2 : Terms in the computation p which cancel. 


commuting s. In the following sections, we will see other expressions of 
this simple idea. 

Algorithm 4.1 can be extended in several different directions. 

1 . The examples above concern vector fields defined on It is possible 

to work out similar algorithms for vector fields defined on more general 
objects, such as Lie groups. This is important for applications in 
robotics and rigid body dynamics. For example, the group G = 50(3) 
is the appropriate configuration space for a rotating rigid body. To be 
more specific, assume the vector fields are of the form 

m=i 

where the are left-invariant vector fields on G and the a ** are smooth 
functions on the group. In this case, the Cayley algebras are generated 
by ordered , labeled trees {16]. Roughly speaking, the trees are ordered 
since the vector fields Y M no longer commute. 

2. The natural action of differential operators on functions turns the ring 

of functions R into a module. In this same way, the trees have a natural 
action on R , as indicated in the Diagram 5. It turns out [34] that this 
gives R the structure of ET/i-bialgebra, as introduced by Nichols [55] 
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and [56]. These types of algebras are closely related to differential 
algebras. 

3. It is a basic fact that the local properties of the nonlinear system ( 1 ) are 
determined by the algebraic properties of the higher order iterated Lie 
brackets; see, for example, [40]. Unfortunately, due to intermediate 
expansion swell, it is often difficult to compute these using current 
computer algebra systems. It turns out that higher order Lie brackets 
not only involve the cancellation of all terms above the first order but 
also the cancellation of some of the first order terms. Algorithm 4.1 
can be used so that the terms arising in these first order cancellations 
need not be computed. 

5 Exponentials, Lie brackets, and nil flows 

Some differential equations have the property that their flows can be inte- 
grated symbolically in closed form. For example, this holds for differential 
equations of the form (1), if E\ and E 2 are homogeneous in the appropriate 
sense and generate a graded, nilpotent Lie algebra. In this section, we give 
an algorithm which, given an appropriate nilpotent Lie algebra, returns vec- 
tor fields on r n with polynomial coefficients which generate the Lie algebra. 
This leads to an interesting class of differential equations whose flows can 
be integrated symbolically in closed form. At the end of the section, we look 
at several applications of this algorithm. 

By the third fundamental theorem of Lie [63], we know that a nilpotent 
Lie algebra arises from some Lie algebra of vector fields. What is not obvious 
is how to construct such vector fields. Nilpotent Lie algebras of vector fields 
have been used as an important tool in partial differential equations by 
Folland and Stein [21], Rothschild and Stein [60], and Rockland [61]; and in 
control theory by Krener [48], Hermes [41], [42], and Crouch [15], [14]. We 
will see below how they are also a useful tool in developing symbolic-numeric 
algorithms to integrate flows. 

Our goal is to describe a natural representation of nilpotent Lie algebras 
on vector fields on Euclidean space with polynomial coefficients. To define 
this represenation, we define a basts of Hall trees on generators E\ and £2 
recursively as follows: 

b basis elements consist of rooted, binary trees, with all nodes, except 
the root, labeled with E\ , £ 2 , £ 3 , • ■ • satisfying 
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(a) all right children are leaves and labeled with either E\ or E 2 

(b) the sequence formed by the labels of the right leaves at increasing 
distance from the root is nonincreasing 

2. the two rooted trees consisting of a root and a single (left) child labeled 
E t) for i — 1,2 are in the basis and of length 1 

3. if we have defined basis elements t t of length 1, . . r — 1 } they are 
simply ordered so that f, < in case the length of t x is less than the 
length of tj. 

4. a tree consisting of a root with a single (left) child labeled E x is in the 
basis provided that E % s left child is in the basis and of lower order. 

To each such tree corresponds an element £\ in the Hall basis [39] of the 
free, nilpotent Lie algebra on two generators. Figure 6 illustrates how the 
basis of Hall trees leads naturally to a representation of the free nilpotent 
Lie algebra generated by E\ and E 2 on the space of vector fields on R ;V 
with polynomial coefficients. See [22] for further details of the following 
algorithm: 

Algorithm 5.1 Fir r > 1 and say that the free , nilpotent Lie algebra of 
two generators of rank r has dimension i V. Let E it . . ., Ejg denote the Hall 
basis. Then the vector fields on R w defined as in Figure 6 satisfy 

1 • the Lie algebra they generate is isomorphic to the free , nilpotent Lie 
algebra on two generators of rank r 

2. any trajectory of the nonlinear system 

x(t) = Ei(x(t)) + u(t) E 2 (x(t)) 9 z(0) = x° 6 R" 

can be computed in closed form in tenns of quadratures of the function 
t^u(t) 

3. there exist constants a\, .... ct/v such that 

i 

( N 

exp(£’ 2 )exp(£' 1 ) = exp 1 £7. 

Vi= 1 

and the a, can he computed by solving a lower triangular linear system. 
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Figure t>: The vector fields arising from the basis of Hall trees. The first 
tree is sent to D \ , while the sum of the trees is sent to D 2 — z\D 3 
+ xix 2 D 5 - ^x^D 6 - \i\x 2 D 7 -iijx^Dg. 


/ 
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As a note, this algorithm was found only after several months of experimen- 
tation with Maple and required the careful study of the free, mlpotent Lie 
algebra on two generators of rank 5, which is 23 dimensional. 

We conclude this section with some remarks describing several applica- 
tions of this algorithm and related work. 

1. Locally approximating a nonlinear system by an explicitly integrable 
mlpotent one yields a number of integration algorithms, which we re- 
fer to as piecewise mlpotent integration algorithms (PWNI). The basic 
idea [23] t [24] is to approximate locally a nonlinear system at a given 
point by an explicitly integrable nilpotent one in which the computa- 
tions can be done symbolically in closed form, and to patch together 
the various mlpotent approximations at nearby points using a stan- 
dard numerical algorithm. Preliminary work indicates that this leads 
to symbolic-numeric algorithms for the path planning problem and ef- 
ficient algorithms to integrate neighborhoods of trajectories around a 
given fixed reference trajectory. 

2. Algorithm 5.1 also provides an efficient means for computing the con- 
catenation of flows. Write x(f) = exp Ex 0 for the flow of the nonlinear 
system 

x(<) = £(*(<)). x(0) = 1 °, (6) 

The Campbell-Baker-FIausdorff formula expresses the concatenation 
of two flows as a single flow: 

exp(<£ 2 ) exp(f£i ) = exp(* E 2 + tE\ + 1/2 t 2 [E 2) E\] 


+ l/12[[£: 2l £7 i], £7i] - 1/12[[£ 2 . £i], £ 2 ] + •••)■ 


Consider the equation 


(7) 


N 

exp(E 2 )exp(E { ) = C ' E <' 

where E t are the vector ^fields produced by the algorithm. Since all 
flows of the vector fields E\ and £2 are explicitly integrable in closed 
form, this reduces the computation of the C{ to the solution of a linear 
system, which is lower triangular. 

3. VVe can also use Algorithm 5.1 to derive a class of numerical integra- 
tors, which are sometimes known as splitting methods. For example, 
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suppose that E\ and E 2 are separately integrable in closed form, but 
E\ -f E'l is not. Then using the algorithm, we can compute constants 
r, such that 


exp(r 7 £ 1 ) • exp(r 6 £ 2 ) ■ exp(r 5 £;i) ■ exp {r 4 E 2 ) 

•exp(r 3 £i) ■ exp(r 2 E 2 ) exp(nEi)- = exp(Ej + £ 2 ) + 0(t 5 ). (8) 

This formula yields a numerical integrator for our original nonlinear 
system (1) (with u = 1). This algorithm is used in accelerator physics 

[65]. 

4. Recently, Strichartz [69] has shown that the solution of the initial value 
problem 

x(t) = E(t,x{t)), x(0) = 1° 

can be written as 

z(t) = exp(G(0) x °i 

with 

T— 1 O 

J I" (^(M 1), £(M2)] ' ' ds < 

and where cr ranges over the symmetric group on r symbols, the inte- 
gration region is a simplex in R r , E(s) denotes the vector field E(s, •), 
and exp is defined as in Equation 6. Formulas of this type date back 
to Chen [13] and appear to be related to Algorithm 5.1. 

5. F. Bergeron, N. Bergeron, and A. M. Garsia have also exploited a rela- 
tion between trees and polynomials in their study of free Lie algebras; 
see [4] and the references cited there. 

6 Taylor series and intrinsic integrators 

Although nonlinear systems often conserve quantities such as energy or an- 
gular momentum, most numerical integrators do not. Similarly, nonlinear 
systems typically evolve on some underlying geometric space, such as a Lie 
group or homogeneous space, but most numerical integrators do not remain 
in such a space. 

Recently, there has been a flurry of activity related to numerical in- 
tegrators preserving the symplectic structure, sparked off by the work of 
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Channel! [10] and Scovel [11], and based upon earlier work of deVogeleaere 
[19]. These types of numerical schemes have found important applications in 
the long term study of orbits in accelerator physics and in other areas [65]. 
The derivation of these integrators typically involves the symbolic computa- 
tion of Taylor series and generating series for the symplectic transformation 
which is the update for the numerical integrator. 

Consider a numerical integrator for a differential equation 

z(t) = E(x(t )), x(0) = p£M 

evolving on a space A/. Call a numerical integrator intrinsic in case x n £ A/ 
implies r n+1 £ Af , for n > 1, where r n is the approximation to the trajectory 
j( 0 at time t n . One means of deriving intrinsic numerical integrators is 
to mimic the derivation of standard numerical integrators, but to impose 
additional constraints on the scheme to satisfy the added condition that the 
points x n remain in the space A/. This typically involves the careful study 
of the Taylor series of the solution. 

This can be done by using the Cayley algebra of trees, as briefly indicated 
in [32], As an illustration of this, we consider intrinsic Runge-Kutta type 
algorithms evolving on a Lie group G, following [16]. Let g denote its Lie 
algebra, and let Y \ , . . Y /v denote a basis of g . We give an algorithm to 
approximate solutions to differential equations evolving on G of the form: 

x(o = £(z(0), x(0) = P eG, 

where 

M=1 

and the = are analytic functions on G. Let exp(t£) ■ x denote the 
solution z(t) at time <. The algorithms depends upon constants c, and c tJ , 

for i = 1, .... £ and j < {. For fixed constants, define the following elements 
of the Lie algebra g 

N 

<i=i 

N 

Ei = ^a^exp (hc2iEi)-z n )Y„eg 

<i = l 
N 

= 5Z a ' 1 ( ex P( /lc 32^2) •exp(hc 31 £ 1 ) • G g ■ ■ ■ 

M=l 
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Figure 7 Trees associated with third order terms in a Taylor series. 

These arise by ‘freezing the coefficients” of E at various points along the 
flow of E . 

Algorithm 6.1 Given an initial point xq on the group , define 
x n + l = exp hckEk ■ • • exp hc } E\z n , 

for n > 0 

Notice that if we assume the exponential exp(/i£ 1 ) maps the Lie algebra 
to the Lie group exactly, then this algorithm is intrinsic. For a group such as 
G = 50(3), there are classical closed form expressions for the exponential 
map and Algorithm 6.1 yields an intrinsic integrator. Notice also that if G 
is the abelian group R N , then the algorithm becomes the classical Runge- 
Kutta algorithm. 

The first step is to derive the equations that the coefficients c, and c i; 
must satisfy in order for the algorithm to yield an rth order numerical in- 
tegrator. This can easily be done using the Cayley algebra of ordered trees 
[16]. The trees are ordered since the vector fields Y M do not commute. 

Assume for the moment that G = and consider the terms in the 
Taylor series 

z(t + h)-x(t) = hz(t) + ^i(t) + + ■ • • 

h 2 

= hE + —DEE 

h 3 f t \ 

+— [DE ■ DE • E + D 2 E{E, £)) + •••. 

Notice that there is a natural correspondence between trees and terms in 
the series. For example, the /i 3 /3! terms are associated with the two labeled 
trees in Figure 7. This observation goes back to at least Cayley [8], [9], 

We now generalize this to a Lie group following [16]. Recall that Dia- 
gram 5 induces an action of trees on the ring of analytic functions on G. 
Using this action, we can now state the 
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Lemma 6.1 Let a denote the tree consisting of a root with a single child 
labeled F. Then for any analytic function f on the group and for sufficiently 
small t, 

f(exp(tF) ■ x ) = exp(fa) • / | x 

Notice that if G is Euclidean space, and if the functions / are the coordinate 
functions x i , . . x n , then this becomes the familar Taylor series. 

Using this lemma, it is now easy to compute the equations that the 
coefficients c x and c x) must satisfy. In spirit, this is similar to Butcher’s use 
of trees to analyze higher order Runge-Kutta algorithms in Euclidean space 

[6], [7]. 
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1 Introduction 


This extended abstract further develops the algorithms in [8] and [9] for 
rewriting expressions involving differential operators. The differential oper- 
ators that we have in mind arise in the local analysis of nonlinear dynamical 
systems. In this work, we extend these algorithms in two different directions: 

1. YVe generalize the algorithms so that they apply to differential opera- 
tors on groups. This generalization is important for applications. For 
example, the nonlinear system describing a robotic joint or a satellite 
evolves on the group G = 50(3) of spatial rotations. The local study 
of such systems requires the computation of expressions consisting of 
differential operators on G. 

2. We develop the data structures and algorithms to compute symboli- 
cally the action of differential operators on functions. Again, this is 
crucial for applications. We illustrate this by deriving conditions for a 
numerical algorithm to remain constrained to a group. In other words, 
if x n +i = T(x n ) is the update rule for a numerical algorithm evolv- 
ing on a group G, we would like to choose T so that x n £ G implies 

■Cn+l € G. 

For a further discussion of applications of these algorithms, see [5] and [6] 
and the references given there. 

Here is the 

Set Up. 

1. Let Ic denote either the real or complex numbers. 

*2. Let G denote a finite dimensional Lie group over k, g denote its Lie 
algebra, and Y\, . . T/v a basis for g of left-invariant vector fields. 

3. Let R = C tro (G) denote the algebra of smooth functions on G taking 
values in k. 

4. Fix M derivations of R of the form 

/ 

N 

F > = T, a j Y »' ( 1 ) 

and let A denote the free associative algebra k<F\, .. ., F\f> of dif- 
ferential operators generated by F\ ) .... Fm, with coefficients from 
k. 
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YVe are concerned with the following 

Problem. Given a differential operator p £ A and a function / £ /?, 
substitute the Equations ( 1 ) and compute p • / using as few operations as 
possible. This problem is interesting since in many cases cancellations take 
place. 

Example 1. Let G — R /V denote the abelian group, 


J dxj 


J = 1 AT, 


the (left invariant) coordinate vector fields, and F\ , F2, F3 three fixed vec- 
tor fields defined in terms of the Y }1 via Equations (1). Then the naive 
substitution of ( 1) and simplification of p • /, where 


P— F3F2F1 — F3F1F2 — F2F1F3 + F1F2F3 £ /l, /£ /i, 

yields 24 A r3 terms, while more specialized algorithms need only compute the 
6/V 3 terms which don't cancel. These types of examples are considered in 
[8] and [ 9 ]. 

Example 2. Consider the local analysis of a nonlinear system of the form 

i (0 = f(z( 0 ), z( 0 ) = z°eG, ( 2 ) 


f = i Z«jFj. 

J=l 

In practice, the u ; are constants, functions of time, or perturbation parame- 
ters. The study of this system typically involves the computations of various 
series in the algebra A of differential operators. For example, the local flow 
of the system is determined by the Taylor series 

exp hF = 1 + hF + f ^F 2 + ^F 3 + . . . € A[[h}}. 

An alternative to computing higher derivatives F k is to choose constants 
c n C{j , i = 1, . . k, j < 2, so that the expression 

exp hckFk • * • exp hc\F\, 
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where 


Fi = X^f* 0 )^ € g 

;i=I 

.V 

F 2 = £ a /i (exp(/tc 2t A) • -r 0 )^ G g 

<i=J 

/V 

F 3 = 5Z « M (exp(/»C 3 2 F 2 ) • exp(/ic 31 F, ) • x 0 )^ G g, 

M= 1 


is equal to exp/iF to order fc. Notice that the loft invariant vector fields F ; 
arise by “freezing the coefficients' 1 of F at various points along its flow. 

Expanding these expressions around the common base point x° £ G 
yields many terms, which must cancel in the end if the algorithm is going 
to approximate the flow of the underlying nonlinear system. The action 
of the differential operators Fj on the coefficient functions a£ must also be 
computed. Notice, that unlike Example 1, the here do not commute. 
This example will be considered in more detail in Section 4. 

The computations in both examples are easily kept track of by using 
finite rooted trees, labeled with the symbols F\ , . . ., Fj^f. It turns out the 
the vector space, with basis the set of such trees, has an algebraic structure 
0 which is crucial to efficiently organizing the computation. The advantage 
of working with the trees B is that many terms which cancel in the end need 
not be computed. See [6] for an expository treatment of this idea. The key 
observation required for this work is that it is possible to define an action of 
the algebra B of finite rooted trees, labeled with F\ 1 . . ., Fm, on the ring of 
functions R which enjoys essentially all the properties of the familiar action 
of the algebra A of differential operators generated by F \ , . . ., F\f on R. It 
turns out that B is a Hopf algebra, just as A is, and that both actions give 
R the structure of what is called an //-module algebra. 

In Section 2, we review the relevant material from algebra. This material 
may be skimmed on a first reading. In Section 3, we define the Hopf algebra 
of Cayley trees and its action on the ring of functions R . In Section 4, we 
continue the discussion of Example 2. 
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2 //-module algebras 

In this section we review the basic facts about bialgebras and //-module 
algebras which will be used in the remainder of this paper. 

In this section, k can be any field of characteristic 0. By an algebra we 
mean a vector space A over the field k with an associative multiplication and 
unit. The multiplication can be represented by a linear map /i : A® k A - A\ 
the unit can be represented by a linear map k — A (the map sending 1 £ k 
to 1 £ /l). The facts that the multiplication is associative, and that 1 £ A 
is a unit, can be expressed by the commutativity of certain diagrams. For 
example, the commutativity of the diagram 

A A ®k A — * A ®k A 

1 l 

A ®k A — - A 

where the upper horizontal arrow is the map fi® /, the left vertical arrow is 
the map / ® /x, and the remaining two arrows are the map /j, expresses the 
associativity of multiplication. 

The dual notion to an algebra is a coalgebra : a vector space C over the 
field k together with a coassociative coproduct A : C — C®kC and a counit 
€ : C — ■ ► k. The fact that A is coassociative and that e is a counit is expressed 
by diagrams which are dual to the diagrams which express the facts that the 
multiplication of an algebra is associative, and that 1 £ A is a unit: they are 
the same diagrams, with the direction of all arrows reversed. For example, 
coassociativitv is expressed by the commutativity of the diagram 

C ®k C ®k C 1 — C ® k C 

r r 

c® k c — c 

where the upper horizontal arrow is the map A ® /, the left vertical arrow 
is the map / ® A. and the remaining two arrows are the map A. Often the 
element A(c) £ C ® k C is written ]C c (i) ® c (2)- 

, < c ) 

A 6ta/ge6ra is a vector space // over k which has both an algebra and 
a coalgebra structure, such that the coalgebra structure maps are algebra 
homomorphisms, or equivalently, the algebra structure maps are coalgebra 
homomorphisms. (This equivalence can be seen by expressing the assertion 
that the coalgebra structure maps are algebra homomorphisms as a set of 
commutative diagrams: this set of diagrams is self-dual.) 

Some examples of bialgebras are the following: 
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1. Let G be a group, and ict kG be the group algebra of G: the vector 
space kG has the elements of G as a basis, with multiplication defined 
by extending the multiplication on G linearly. The coproduct and 
counit of kG are defined by 


Mg) = <j®g \ 
<(9) = 1 j 


g G G. 


2. Let G be an affine algebraic group, and let k[G] be the algebra of 
representative functions on G. The algebra structure of k[G ] is the 
usual algebra structure of functions with point-wise multiplications. 
The coproduct arises from the group multiplication G x G — ■ G, which 
induces the map k[G] — k[G x G] ~ k[G] k[G\. The counit arises 
from the map {e} — * G % where {e} is the single-element group. 

3. Let L be a Lie algebra over k , and let U(L) be the universal enveloping 
algebra of L. The coproduct and counit of U(L) are defined by 


A(x) = 1 ® x + x ® 1 

t(x) = 0 


xGi, 


and extended to all of U(L) using the fact that A and e are algebra 
homomorphisms. 

Usually, in studying bialgebras, an additional condition is imposed which 
is analogous to the assertion that a semigroup is a group. Such bialgebras are 
called Hopf algebras. The bialgebras which we consider in this paper (such 
as the universal enveloping algebra of a Lie algebra) satisfy this condition 
automatically. 

A coalgebra is said to be cocommutative if it satisfies A = To A, where 
T is the map T : C C —> C ®* C defined by T{x ® y) = y ® x. Note that 
the bialgebras in Examples 1 and 3 are cocommutative. 

A vector space V over k is said to be graded if there is a sequence of 
subspaces Vo, V\, ... such that / 


v = ® v n . 

n= 0 


A graded vector space V is said to be connected if Vo = k. 
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Let H be a bialgebra. A II -module algebra is an algebra ft which is an 
//-module such that the action satisfies 

/> • (/<7) = X^(/»( 1 ) /)(/>( 2 ) • !/), for all hell, f,g£R. 

<M 

Remark 2.1 If g (E // satisfies A(g) ~ g®g and R is an //-module algebra, 
then g acts as an endomorphism of R\ if x £ // satisfies A(x) = l <g> x + x 0 1 
and R is an //-module algebra, then x acts as a derivation of R. 

3 //-module algebras and Cayley trees 

In this section we describe a bialgebra structure on the vector space with 
basis all equivalence classes of rooted trees. The relation between trees 
and differential operators goes back at least as far as Cayley [3] and [-1] . 
Important use of this relation has been made by Butcher in his work on 
higher order Runge-Kutta algorithms [1] and [2], In this section and the 
next, we follow the treatment in [8] and [9]. By a tree we mean a nonempty- 
finite rooted tree, and by a forest we mean a finite family of finite rooted 
trees, possibly empty. 

Suppose { F \ , ... F\f } is a set of formal symbols (which later will be the 
names of differential operators). By a labeled tree we mean a tree for which 
we have assigned an element of { Fi , .... F M } to each node, other than 
the root, of the tree. We say that a tree is ordered in case there is a partial 
ordering on the nodes such that the children of each node are non-decreasing 
with respect to the ordering. 

We now describe the bialgebra structure on spaces of trees. Let 

HT(F u .... F M )} 

denote the vector space which has as basis all equivalence classes of labeled, 
ordered trees. The vector space k{T{F \ , ..., F w )} is graded, with the grad- 
ing given as follows: if the tree t has n + 1 nodes, then 

t G • • • , F \ * )}„. 

We now define the multiplication on k{T(F u ..., Fa/)}. Since the set of 
labeled, ordered trees form a basis for *{T(F, , ..., Fa/)}, it is sufficient to 
describe the product of two such trees. Suppose that f, and t 2 are labeled, 
ordered trees. Let si , . . . , s T be the children of the root of fj . If t 2 has n + I 
nodes (counting the root), there are (n -(- l) r ways to attach the r subtrees 
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of t\ which have sj , . . . , s r as roots to the labeled tree <2 by making each 
.s, the child of some node of < 2 , keeping all the original labels. Order the 
nodes in the product so that the nodes which originally belonged to each tree 
retain the same relative order to each other, but all the nodes that orginally 
belonged to <1 are greater in the ordering than the nodes that originally 
belonged to <2 The product t]t 2 is defined to be the sum of these (n + l) r 
labeled trees. It can be shown that this product is associative, and that the 
trivial labeled tree consisting only of the (unlabeled) root is a right and left 
unit for this product. For details, see [7]. 

We now define t he comultiplication on k{T(F x , . . . , Fm)}. If t is a tree 
whose root has children s,,...,s r , the coproduct A (t) is the sum of the 2 r 
terms 1 1 Qf 2 , where the children of the root of t, and the children of the root 
of t 2 range over all 2 r possible partitions of the children of the root oft into 
two subsets. I he labels remain the same, and the ordering is handled in the 
same wav as in the product. The map e which sends the trivial labeled tree 
to 1 and every other tree to 0 is a counit for this coproduct. In [7], it is 
shown that these algebra and coalgebra structures are compatible, proving 
the 

Theorem 3.1 /Tie space k{T(F\, . . ., F\f)} is a graded connected cocom- 
mutative bialgebra. 

We call this algebra the algebra of Cayley trees. 

We now define an action of the algebra of Caylev trees 

Z? = k{T(F u .... F m )) 

on the ring R. making R a S-module algebra, which captures the action of 
trees as higher derivations. The action is defined using the map 

V : k{T(Fi, • * , Fm)} — End*. /?, 

as follows: 

1. Given a labeled, ordered tree t with m + 1 nodes, assign the root 
the number 0 and assign the remaining nodes the numbers 1, . . . , m. 
We identify the node with/the number assigned to it. To the node k 
associate the summation index /i*. Denote (/ii , . . . , ^ m ) by p. 

2. For the labeled tree f, let k be a node off, labeled with F^ k if k > 0, 

and let / V be the children of k . Define 

R(k;fi) = Ym ■ * ■ Y ^ t , , if k > 0 is not the root; 

— ■ • * Yu ., , if k = 0 is the root. 
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Note that if* > 0, then R(k\p) € R. 


3. Define 

i v 

^(0 = H /?(m;/i) • • • ft(l;/i)c(0;/i). 

Mm=l 

4. Extend xp to all of lc{T(F \ , . . . , f\/)} by linearity. 

It is straightforward to check that this action of D on R makes R into a 
5-module algebra. 

We summarize with the following theorem. 

Theorem 3.2 Let G denote a finite dimensional Lie group and R the alge- 
bra of smooth functions on G, as detailed in the Set Up. Let 0 denote the 
algebra of Cayley trees k{T(F \ , .... F\i)}. Then R is a O-module algebra 
with respect to the action defined by xp. 

Remark 3.1 The standard action of the algebra /I of differential operators 
generated by F i , . . ., F\j on the algebra of smooth functions R gives R the 
structure of a /1-module algebra. It is easy to relate these two //-module 
algebra structures on R and this observation is the basis for our algorithms. 
Let 

4> : A — . B 

denote the map sending the generator Fj of the algebra A to the tree con- 
sisting of two nodes: the root and a single child labeled Fj. Extend 0 to be 
an algebra homomorphism. Let \ denote the map 

/! — End* R 

defined by using the substitution (1) and simplifying to obtain an endomor- 
phim of R. 


Theorem 3.3 (i) The maps y, 0 and \p are related by y = xpoo. (u) Fix 

a function f £ R and a differential operator p £ A. Then 

/ 

P * / = <P(P) * f 

Here the action on the left views R as an A-module algebra , while the action 
on the right views R as B -module algebra. 
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The first assertion is proved in [9] and the second assertion follows from the 
first assertion and the definitions. 

Using this theorem, it is easy to give an algorithm to solve the Problem 
posed in Section I. We defer to later paper a complete analysis of the 
complexity of the algorithm and simply remark here that in many examples 
the algorithm results in a savings which is exponential in the degree of the 
differential operator. 

Algorithm. Given a smooth function f £ R and a differential operator 
p £ /l, compute the function p • / via <p(p) . /. 

4 Applications 

We use the notation of the Set Up from Section 1 . Let exp(/iF)x denote the 
resulting of flowing for time h along the trajectory of the nonlinear system 
(2) through the initial point x° £ G. We require a theorem concerned with 
the explicit computation of terms in the Taylor series expansion of a solution 
of (2). This is one of the main applications of the symbolic calculus described 
in the sections above. 

This theorem is most easily stated if we introduce two additional oper- 
ations on the algebra of Cayley trees B. Given a,0 <= B, define the meld 
product 0 0 or to be the labeled, ordered tree obtained by identifying the 
roots of the two trees. The meld product is then extended to all of B by 
linearity. Given a derivation F 6 Der(fl), let 0 be the tree <p(F) and let 
or e B. Recall 0 is a tree consisting of a root and a node labeled F. We 
define the composition product Boa to be the tree formed by attaching the 
subtrees whose roots are the children of the root of a to the node labeled F 
of the tree 0. If a € B is a tree, define the exponential and Meld-exponential 
of a tree by the formal power series 

exp(Aar) = 1 + ha + L-r* 2 + —a 3 + ■■■ 

2 . 3 ! 

Mexp( ha) = 1 + /iar + — a 0 a + — a 0 QrQ a ^_... 

Theorem 4.1 (i) Assume f <= R and F £ Der(ff). If f is analytic near x, 
then for sufficiently small h, 

f (exp(hF)x) = exp {h<f>{F)) ■ f\ x . 
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(») Let F = Em=i a"(exp(/iG)x u )y )11 where G £ Der(ft), and x° £ G. 
Then 


F ■ f = (<p(F) o Mexp(hG)) • /. 


Using this theorem, it is easy to analyze the numerical algorithm de- 
scribed in Example 2 of Section 1. For typographical reasons, we use the 
following one dimensional notation for trees 1 : the tree consisting of a root 
and a single child labeled F\ is denoted /[Fj]; the tree consisting of a root 
and two children labeled F\ and F2 is denoted /[Fj,/^]; the tree consisting 
of a root, with a single child labeled F\ , which itself has two children labeled 
F2 and F3 is denoted I[F\ [F2 , Fi]], etc. Note that the labels need not be 
distinct, but their order is important. 

Consider the expression 


exp /1C3 /' '3 exp I1C2 /'2 exp hc\ F\ 


computed to order h 3 . Let p G A denote the resulting expression. The 
image <p(p) G B contains the following terms: 


h 3 


c 2 c 2 l 

2 ! 


HF[F, F}] + 


/. 3 


—■nnr, fii + ^ 


3C32 


2 ! 


/[F[F,F]] 


+ /» 3 C 3 C 3 |C 3 2 /[F[F, F]j. 


Our goal is to choose the constants c, and c X] so that that 


exp /iF = p + 0 (h 4 ). 


One of the third order term arising from <p(exphF) is ^-/[F, [F, F]]. Setting 
the coefficients of these trees equal to each other yields the constraint: 


cjcji 

2 ! 


c 3 c ai 

2 ! 


+ 


<*<4 . _ 1 

— jjj •" C 3C31C32 — 


Other constraints arise from the other trees. We have coded this algorithm 
in Maple, Mathematica, and Srtobol and are currently experimenting with 
it. 

The notation is due to Peter Olver, as is some of the Mathematica code used to 
generate these examples. 
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Bialgebra deformations and algebras of trees* 

Robert Grossman* and David Radford 


Abstract 

Let A denote a bialgebra over a field k and let A t = A[[£]] denote 
the ring of formal power series with coefficients in A. Assume that A is 
also isomorphic to a free, associative algebra over k. We give a simple 
construction which makes A t a bialgebra deformation of A. In typical 
applications, A t is neither commutative nor cocommutative. In the 
terminology of [1], A t is a quantum group. This construction yields 
quantum groups associated with families of trees. 

1 Introduction 

Let A denote a bialgebra over a field k and let At = A[[f]] denote the ring of 
formal power series with coefficients in A. Assume that A is also isomorphic 
to a free, associative algebra over k. We give a simple construction which 
makes A t a bialgebra deformation of A . In typical applications, At is neither 
commutative nor cocommutative. In the terminology of [1], At is a quantum 
group. This is an extended abstract. The final detailed version of this 
paper [6] containing proofs and additional examples has been submitted for 
publication elsewhere. 

An interesting class of examples is obtained by taking the bialgebra A 
to be the Hopf algebra associated with certain families of trees as in [4]. 
In fact, these examples are closely related to each other and to algorithms 
pertaining to differential operators [5]. 

Formal deformations of bialgebras and quantum groups has also been 
studied from related viewpoints by Drinfeld [1], Gerstenhaber [2], and Ger- 
stenhaber and Schack [3] 

*This is an extended abstract. The final detailed version of this paper has been sub- 
mitted for publication elsewhere. 

f This research is supported in part by NASA grant NAG2-513, NSF Grant DMS- 
8904740, and by the Laboratory for Advanced Computing at the University of Illinois at 
Chicago. 
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2 Power series and the completed tensor product 

Suppose that k is a field and A and B are it-algebras. We let At = /![[£]] 
denote the ring of formal power series over A with its usual multiplication. 
Let / : A — ► B t be a fc-linear map and write f(a) = c/(a, n )^ n for 

a G A. Define a /t-linear map / : At — * Bt by 

^ OO oo 

/( a ) = c /( a iJ)) tn > for a = ^ a n t n G A t . 

n=0 «-f j=n n=0 

Observe that we can define a cateogy {Algk) ti whose objects are A t and 
whose morphisms are fc-iinear maps / : A t — * B t satisfying / = / 1^, where 
/\a is the restriction of / to A. 

We define the completed tensor product A t ^^ t B t of A t and B t over k t in 
this category to be (A®k B)t • For a = ££°=o a n t n G A t and b = b n t n G 

Bt we let 

oo 

a<g>b=^T( ^2 ai ®bj)t n . 
n= 0 i+jssn 

Suppose that / : A t — *- A[ and g : Bt — ► B[ are morphisms. We define 

a morphism of completed tensor products f(&g : At®k t Bt *■ A' t ®k,B' t by 

setting (/® 3 )U® k a(a ® b) = f(a)®g(b) for a G A and 6 G B. The usual 
formalism for the linear tensor product of maps translates to 

(f®g){a®b) = /(a)®</(b) 

for a G At and b G B t in this category. Note that if / and g are also algebra 
maps, then the morphism f®g : At®k t Bt — * A' t ®k,B' t is an algebra map. 

Now let C be any coalgebra over k. A sequence of elements cq, c\ , C 2 , ... G 
C is called a sequence of divided powers if 

A(c n ) = "^2 c « ® C J anc * c(c n ) = ^o,n for all n > 0. 

:'+j=n 

Probably the most basic example of a sequence of divided powers arises 
from a primitive element in a bialgebra. Suppose that A is a bialgebra over 
k and that £ G A is primitive. Since A is multiplicative, we calculate by the 
binomial theorem that A(£ n ) = (A(£)) n = (t?®l+l®^) n = 
for n > 0. Thus co, ci,C 2 , ... is a sequence of divided powers, where c n = ^ 
for n > 0. 
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3 Deformations of certain enveloping algebras 

The notions of algebra, coalgebra, bialgebra and Hopf algebra in the category 
{Algk)t are the same as those in the category of vector spaces over the field 
k except, of course, the structure maps are required to be morphisms. Let 

(A, m, 77) be an algebra over k, where m : A ® A * A is multiplication and 

77 : k * A defines the unity of A. Then ( A m, fj) is an algebra in (Alg^t. 

It is easy to see that rh(a®b) = ab for a, b £ At . A morphism / : At — *■ B t 
is a morphism of algebras if and only /\a : A — * B t is a map of t-algebras. 

The proof of the proposition below is really a matter of unravelling def- 
initions. 

Proposition 1 Suppose that A is an algebra over a field k with a k-coalgebra 
structure (A, A,c). Then (/l*,A,c) is a coalgebra in (Alg^t- Furthermore 

^ 00 00 

A(») = atld ? ( a ) = T, f ( a n) <n 

n-i) n=0 

for a = a n t n G A t . 

Suppose that ( A t , A,e) is a coalgebra in (Algk)t- We say that K G A t is 
grouplike if 

A(I<) = K%K and e(I<) = 1. 

We say that £ € A t is nearly primitive if 

A(£) = t®I< + H®£ 

for some grouplike elements K, H £ A t . If K = H = 1 then £ is said to be 
primitive. 

For an algebra A over a field k of characteristic 0 we let exp(at) = 
H^Lo{^r)l n G A t . The following corollary gives the relationship between 
sequences of divided powers and grouplike elements. 

Corollary 1 Suppose that A is an algebra over a field k which has a k - 
coagebra structure (T,A,e). Let (At, A,?) be the resulting coalgebra in 
( Algk)t • Then: 

(a) Let K = Y^Lo a nt n G A t . Then K is grouplike if and only if do, a .\ , a 2 , ... 
is a sequence of divided powers in A. 

(b) Suppose that the characteristic of k is 0 . If a £ A is primitive, then 
K = exp(at) is a grouplike element of At, 
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Now we construct deformations of enveloping algebras over a field of 
characteristic 0 which are free as associative algebras on a space of primitives. 

Theorem 1 Suppose that V is a vector space over a field k of characteristic 
0. Turn the tensor algebra of V into a bialgebra (j 4, m, r\ ) A, e) by 

defining A(£) =£®l+l®£ and e(£) — 0 for £ € V . Let p,q G V and write 
V as a direct sum of subspaces V — P © P f , where P = span(p, q). Then 
there is bialgebra deformation (A t , m, f;, A, e) of (A, m, rj, A, e) such that 

a) A{£) = £®l + l®£ for £ G P, 

b) I\ — exp(pt) and H ~ exp(qt) are grouplike elements of (At, A, ?), 
and 

c) A(£) = £®K + H®£ for £ £ P* . 

We comment that I\ — exp(tp) = 1 when p = 0. If p ^ 0 and <7 = 0, 
for example, then the deformation of the theorem is not cocommutative. If 
dim(V) > I then the free algebra A is not commutative. In this case the 
deformation of the theorem is not commutative. 

4 Deformations of bialgebras trees 

In this section we give an example from [4] involving bialgebras of trees to 
which Theorem 1 applies. Let k be a field of characteristic 0. Let T be the 
set of finite rooted trees, and let k{T } be the fc-vector space which has T 
as a basis. 

We first define a coalgebra structure on k{T}. If t £ T is a tree whose 
root has children s\ , . . . , s r , the coproduct A (t) is the sum of the 2 r terms 
® *2i where the children of the root of t\ and the children of the root of 
range over all T possible partitions of the children of the root of t into two 
subsets. The map e which sends the trivial tree to 1 and every other tree 
to 0 is a counit for this coproduct. It is easy to see that comultiplication is 
cocommutative. 

We now define an algebra structure on fc{T}. Suppose that fj, £ T 
are trees. Let s i, . s r be the children of the root of t\. If t<i has n + 1 
nodes (counting the root), there are (n + l) r ways to attach the r subtrees of 
1 1 which have . . . , s r as roots to the tree *2 by making each s, the child of 
some node of t 2 . The product t\t 2 is defined to be the sum of these (n-b l) r 
trees. It turns out that this product is associative, and that the trivial tree 
e consisting only of the root is a right and left unit. It can also be shown 
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that the maps defining the coalgebra structure are algebra homomorphisms, 
so that A = k{T} is a bialgebra. 

For technical reasons, we require that nodes of the tree be ordered. We 
say that a rooted finite tree is ordered in case there is a partial ordering 
on the nodes such that the children of each node are non-decreasing with 
respect to the ordering. To define the product of two ordered trees t\ and 
t 2 j compute the product of the underlying trees, and order the nodes in the 
product so that the nodes which originally belonged to each tree retain the 
same relative order to each other, but all the nodes that orginally belonged 
to are greater in the ordering than the nodes that originally belonged to 
^2 • 

Let X be the set of of trees whose root has one child. Then A =. 
k< I(X)> as an algebra [4]. Applying Theorem 1 yields a bialgebra de- 
formation A t of A which is neither commutative nor cocommutative. 
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